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NOMENCLATURE 

a  radius  of  the  circular  cylinder 

a.|j  added  mass  or  moment  coefficient 

d  gap  between  the  cylinder  and  the  bottom  of  the  sea 

0  volume  of  the  cylinder 

g  gravitational  acceleration 

h  depth  of  water 

H|^  J|^  +  1  Y|^  Hankel  function 


I|^  Modified  Bessel  function,  first  kind  of  order  k 

Bessel  function  first  kind  of  order  k 

K  Modified  Bessel  function,  second  kind  of  order  k 
n 

mo,mq  as  defined  by  the equation  2.22,  2.23 
n  normal  vector 

?  position  vector 

t  draft  of  the  cylinder 

V  velocity 

x,y,z  cartesian  coordinate  system 
Y|(  Bessel  function,  second  kind  of  order  k 
Z|^  function  defined  by  equation  2.18,  2.19 
Kronecker  delta 

0  angular  position  In  cylindrical  coordinates 
p  fluid  density 

velocity  potentials 

V  2,,. 


radian  frequency 


I  OUxr 

TAil 

I  ^lajaounoed 

I  Avail t 
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HYDROOYHAMIC  COEFFICIENTS  FOR  VERTICAL  CIRCULAR  CYLINDERS 

AT  FINITE  DEPTH 


Vertical  circular  cylinders  are  used  In  many  oceanographic  applications 
such  as  buoys,  drilling  rigs  and  instrumentation  platform  for  their  simplicity 
In  construction.  The  available  data  on  their  hydrodynamic  coefficients  Is 
limited  and  the  present  numerical  procedures  are  based  on  finite  element 
solutions  or  the  numerical  solutions  of  Integral  equations.  The  present 
method  Is  relatively  simple  to  formulate  and  the  solution  requires  a  very 
short  computer  time.  The  hydrodynamic  coefficients  such  as  added  mass, 
damping  coefficients  for  heave,  sway  and  pitch  motions  are  formulated  and 
the  results  for  different  depth  to  radius  and  draft  to  radius  values  are 
presented  In  graphical  form. 

Havelock  (1955)  theoretically  determined  the  added  mass  and 
damping  coefficients  for  a  sphere.  Kim  (1965)  studied  the  hydrodynamic 
coefficients  for  ellpsoldal  bodies  oscillating  at  the  free  surface.  Shen 
Wang  (1966)  calculated  the  added  mass  and  damping  coefficients  of  sphere 
In  Infinite  and  finite  depth  of  water.  Garrison  (1975)  gave  the  general 
formulation  of  these  coefficients  for  arbitrary  forms  In  terms  of 
distributed  singularities  and  the  numerical  results  for  a  vertical  circular 


cylinder  in  infinite  and  finite  depth  of  water.  Bai  and  Yeung  (1974) 
calculated  added  mass  coefficients  for  horizontal  and  vertical  cylinders. 
Bai  (1976)  gave  the  added  mass  and  damping  coefficients  for  axisymmetric 
ocean  platforms  Kritis  (1979)  had  applied  the  hybrid  integral  method 
of  Yeung  to  axisymmetric  bodies  and  gave  numerical  results  for  a  circular 
cylinder. 

The  various  methods  developed  for  the  solution  of  three  dimensional 
axi symmetrical  bodies  can  be  summarized  as  follows.  In  the  first  group 
of  methods  Sources  and  Multipoles  are  distributed  inside  the  body  and  their 
strength  is  calculated  to  satisfy  impermeable  boundary  conditions  of  the 
body.  The  second  set  of  solutions  distributes  the  singularities  at  the 
surface  of  the  body  and  an  integral  equation  is  used  through  the  use 
of  Green's  theorem.  The  solution  of  the  integral  equation  gives  the 
strength  of  the  singularities.  Thirdlyf  the  finite  element  formulation  is 
used  to  find  the  velocity  potential  at  specified  node  points.  Possible 
combinations  such  as  the  Hybrid  method  referred  above  also  exist,  combining 
the  above  solutions  and  reducing  the  computational  effort.  The  present 
formulation  follows  the  general  procedure  outlined  by  Garrett  (1970)  who 
studied  the  scattering  of  waves  at  the  presence  of  circular  docks. 

Although  it  is  of  major  concern  to  Naval  Architects  and  ocean 
engineers  very  few  data  exist  on  the  hydrodynamic  coefficients  of  circular 
cylinders.  Serving  to  this  aim  graphical  results  covering  a  large  range 
of  parameters  and  sutnnary  formulas  for  computer  applications  are  presented 
in  this  paper. 

2.  FORMULATION  AND  SOLUTION  OF  THE  PROBLEM 

The  coordinate  system  Oxyz  is  shown  in  Figure  1.  The  origin  is 
at  the  bottom  and  ^  is  positive  upwards.  The  region  O^z^h  is  assumed 


to  be  filled  by  an  Incompressible  fluid  of  density  p.  The  undisturbed 
free  surface  Is  at  z  -  h.  The  radius  of  the  cylinder  Is  a  and  draft 
T  =  h  -  di.where  d  Is  the  gap  between  the  cylinder  and  the  bottom.  The 
standard  small  motion  assumptions  are  made  and  the  motion  Is  periodic 
with  frequency  0).  An  irrotatlonal  flow  Is  assumed  to  exist  given  by 
$  (r,0,z;t)  =  Re{${r,0,z)  e“’“S  (2.1) 

where  r,0,z  are  cylindrical  coordinates  and  0  =  0  corresponds  to  the 
positive  X  axis.  ^  {r,0,z)  Is  a  complex  spatial  velocity  potential  which 


satisfies: 


ir  .11 

r  ar 


(2.2) 


In  the  fluid  region. 

The  boundary  conditions  for  different  motions  of  the  cylinder  are 


given  below. 


At  the  free  surface 


r  ®  =  o  atz  =  h 


*  0  at  z  »  0 
-n  “  +  ^‘(^  X  fi) 


(2.3) 

(2.4) 

(2.5) 


on  the  body  surface  where  n  Is  the  normal  to  the  surface.  For  heave  motion 


on  z  =  d 


(2.6) 


~  =  o  on  r  *  a  for  d<  z<  h 
3r  —  — 


(2.7) 


where  Is  the  velocity  of  the  cylinder  due  to  heaving  motion  alone- 
For  sway  motion 


?z‘  “  0  on  z  =  d 


(2.8) 

(2.9) 


I"  *  Vs  COS0  at  r  *  a  for  d<^  z^  h  (2. 

where  Vg  Is  the  velocity  of  the  cylinder  due  to  sway  motion  alone. 
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f 


For  pitching  motion 

*  ft  r  cose  (2.10) 

1“  =  ft  (h-z)  cose  (2.11) 

where  ft  Is  the  angular  speed  of  pitching  motion. 

Following  Gerrett's  method  the  fluid  Is  divided  Into  two  parts  namely 
the  Interior  part  ABCD  and  the  exterior  part  ACEF  as  In  Figure  1.  The 
appropriate  solution  of  velocity  potential  in  each  region  is  found 
and  the  solutions  are  matched  at  the  boundary  so  as  to  have  the  continuity 
In  ^  and  its  first  derivatives  are  satisfied.  The  velocity  potential  In 
the  interior  domain  is  expressed  as 

k  L 

<l>  =  [♦p  (r»z)  (r.z)]  cosKe  for  o<r<a  o<z<  d  (2.12) 

where  k  =  o  refers  to  heave  motion  and  k  =  1  refers  to  sway  and  pitch  motion, 
k  k 

<frp  and  <j>  are  the  particular  and  homogeneous  solutions  in  the  Interior 
region. 

The  particular  solutions  for  different  motions  are  given  as 

*p  *  2^  (z^  -  )  for  heave  (*)  (2.13) 

t  p  (r,2)  »  0  for  sway  (2.14) 

^  p  “  23T  (2.15) 

The  particular  solutions  satisfy  the  respective  kinematic  conditions  at 
the  bottom  of  the  cylinder  and  the  bottom  of  the  fluid  for  r^a  given 

k. 

by  the  equation  (2.4)  through  (2.11).  D  Is  a  constant  of  dimension 
[L  /T]  and  Is  chosen  to  fit  the  particular  motion  of  the  cylinder. 

In  the  exterior  region  the  velocity  potential  Is  given  In  terms  of 
eigen  expansion. 

(*)  A  particular  solution  for  heave  motion  Is  suggested  by  Professor  J.N. 

Newman  of  M.I.T. 
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♦  =  ^  e 


(2.16) 


“k  ^0  ^  K  \  K  ^q 

q*l  *< 

where  Hj^  Is  the  Hankel  function  of  the  first  kind  of  order  K  ®  0,1 

and  K|(  Is  the  modified  Bessel  function  of  the  second  kind  of  order  K  =  0,1. 
k  k 

B*:  and  B^  are  complex  unknowns.  The  orthonormal  Z,.  functions  In  the  Interval 
0  q  1 

o<z<h  are  defined  as 


Zq  (z)  =  No  '  cosh  (m^z) 
Zq  (z)  *  N^^/^  cos  im^z) 


(2.18) 

(2.19) 


No  =  7  [1  + 


sinh  (2  mr.hl 


where  mo  and  mg  are  the  solution  of  the  equations 

J 

mo  tanh  (m^h)  =  g 

2 

fflq  tan  (w^h)  *  -  ^ 


(2.20) 


(2.21) 


(2.22) 

(2.23) 


To  satisfy  the  continuity  of  velocity  potential  at  r  =  a  for 


o<z<d  homogeneous  solution  defined  above  Is  used  as  follows 
♦p  (a,z)  +  (a,z)  =  ^  (a,z) 


J_  [♦!!  (a.z)  +  (a,z)]  =  (a,z) 

3r  P  ® 


(2.24) 

(2.25 


for  o<z<d. 


The  value  of  homogeneous  potential  at  r  -  a  can  be  expanded  as 


♦  h  (a»z)  =  ^  +  An  cos  (^^) 


(2.26) 


For  continuity  of  potential  function  coefficients  a|['s  are  obtained 
using  (2.24)  as: 


^n  “  I  T  Ue  ■  ♦S  («»2)]  cos  (^)  dz  (2.27) 
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The  particular  solution  Is  a  known  function  therefore  can 


be  expressed  as 


*5  •  I  o'  <*•')  (t)  -  % 


(2.28) 


where 


k  9  A  it  k=0.  , 

“  n  *  J  ^  *0  ^a.?)  cos  (nirz/d)  dz 

0 


The  Interior  homogeneous  solution  can  be  written  now  as 


♦  5  C'-^)  ■  i  (J)'  *  j,  <  ^^^0,  cos  (SJi)  (2.30) 


in  r<a,  o<z<d 


Multiplying  the  equation  (2.25)  by  Zq  (z)  and  considering  that 
1^  (a,e,z)  =  V  cose  for  d^z^h  and  integrating  between  limits  In  o  and  h 


one  obtains 


B  J  •  (h  m^)  H,,’  y  l  *  Zo  »  5  (2.31.a) 


d  ,  k 


(mqh)  .  (mqa)  f  *  Zq  U)  dz  +  a  ^ 


forq  =  1,2,.. 


wherf 


“  0  Z^,  (z)  dz  +  pj-  (z)  dz 


Bq  =  /  ^^ar^^  Zq  (z)  dz  +  pj-  /  ^  ^  (z)  dz 


(2.31-b) 


(2.32) 


(2.33) 


Where  V  =  0  for  heave,  V  »  Vg  for  sway  and  V  *  n  (h-z)  for  pitch 


motion. 


Inserting  the  values  of  J  from  (2.17)  Into  equation  (2.27)  and 
4i|5  from  (2.28)  into  (2.31-a)  and  (2.31-b)  we  obtain 
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00 

K  ?  If  n,  (nira/d)  for  n  =  0,1,2,.. 

bJ  *  -C  ^0  •  ^00  ^  n=1  ^  T  Ik  (nira/d)  *  ^on>  &k  (2.35-a) 


Cm^h)  (m^a)  (m^h)  (m^^a) 


c  iv*^  I  ?  (nifa/d) 

.  ^?a  ^0  *  %  ^  ^  T  Ik  {n7a/d) 

Bq»  - -  + 


jJL 


(niqh)  K|^  (m^a) 


(2.35-b) 


(niqb)  K,^’(m^^a) 

The  equations  (2.34),  and  (2.35-a,b}  form  a  coupled  system  of  equations 

where  and  b{$  are  the  unknown^.  Substituting  from  (2.35-a,b)  Into 
n  q  ’  q 

(2.34),  a  set  of  linear  equations  for  are  obtained. 

Thus:  „  k  t  k  .  . 

^  aJ  =  -  hj  for  n  =  0,  1,  2,  ...  (2.36) 

j=o  " 


where 


k  -  -."kK^)  W  hi  .  .  MV’ 

no  (m^a)  m^h  ij,,  K^'ft^a) 


.k  .  ,"k<%»>  “o  •  ^on  .  " 

"n  m^h 


q=l 


3k  E 
g  ♦  an 

mqh 


)  - 


^  (2.37) 

niqh  J  nj  ' 


(2.38) 


where  kronecker  delta  , 
nj 


uj  =  ^  l’  (jTra/d)/l,^(jwa/d), 


for  j  =  1,2,3,. 


,,k  _  kd 
^0  "  4?  ’ 

_  2 


"qn  "  J 


/  Z_(z)cos  (mrz/d)  dz 


(2.39 

(2.40 

(2.41) 


3.  HYDRODYNAMIC  FORCES  AND  MOMENTS  IN  TERMS  OF  THE  VELOCITY  POTENTIAL. 
The  forces  and  moments  are  defined  by  the  Integrals  taken  over  the 
body  surface  as  follows: 
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(3.1) 


//♦•nds 

=  -  p  ^  //  ♦  •  (r  X  ds 


(3.2) 


The  added  mass  and  damping  coefficients  are  calculated  in  the  following 


manner.  aii  hi* 

'TO' ^  ”  P D  ^Mr.e.z)  ds 

(r.G.2)-(r  n)  ds 


(3.3) 


(3.4) 


®66  ^  .  ^66  ,,.l  /  \  ■^\  j 


(3.5) 


(3.6) 


^Mr.e.z)  n^  ds  (3.7) 

where  D  is  the  volume  of  the  cylinder,  a-j^  and  b^-j  are  sway  added 
mass  and  damping  coefficients  and  a-jg  and  b^g  are  sway  induced  pitch 
added  moment  and  damping  coefficient,  a22  and  b22  are  the  heave  added 
mass  and  damping  coefficients,  ag^  and  bg-j  pitch  induced  added  mass 
and  damping  coefficients.  Here  in  (3.3)  and  (3.4)  is  the  potential 
for  sway  motions,  and  4'^  in  (3.6)  and  (3.7)  is  the  potential  for  pitch 
motion. 

For  reasons  of  symnetry 
*16  “  *61 
^16  “  *’61 
*44  *  *66 
*^44  "  ^66 
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The  heave  added  mass  and  damping  coefficients  in  particular  are  given 


Using  4)®  and  41^  form  quation  (2.13)  and  (2.30)  into  (3.8)  and  integrating 

a  K  9  -  IP  f  1  Ii  (nira/d) 

f22  .  .  '’22  _  d  .  1  1  ,aN^  .  1  .0  .  2  ,dx  ^'^=11.  1  -  y  (3.9) 

Similar  calculations  are  made  for  other  added  mass  and  damping 
coefficients  and  these  formulas  are  given  below. 


in  .  ,  in  ,  ,  „i  "i(V>  tsh(%h 

pD  pD  0  N  l/'i  m  (h-d 

0  0 


H^(mQa)  [sh(mjjh)-shmjjd]  ^  gl  (tiLa)[sin(nLh)-sin(rtLd)] 

0  N^l/a  n,„(h-d) - “I  -  r^,  1/2  (h-d) 

(3.10) 


<?  Ct„-N-'^2„^(V)-  bJ  H-f  t,  K,  (h^a)  ^’]  + 


1  ,ax^(l-z-(l-)^)  +  i  (t)  ^  A  ^2  , 

-[ff  (a-)  8-  0  .  n-o  n  %  ]> 


S  /  -inH 


'  SB-a=  ,„^a)  bI*J,  * 

.rl  (1)  a’  *  1  I  a'  n 

LQ  ^  XT  ^  r>  r  J-r 


‘■8  'd'  0  n  ^n  I-j  (rnra/d) 


where 


+•  =  snA'i  sh(mnd)  j.  ch(m{^)-ch(m  4)  , 

^0  - (H^-i - > 

•d  =  <n-J)  -  costoil^^soUsi,!), 


(3.12) 


(3.11) 


<  "J  ",  (V'  -S-i'T'-  *  1  »’<. 


sin  -  sin  "'l^ 


rn^Th-d, 


(3.13) 
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4.  LIMITING  VALUES  FOR  HEAVE  ADDED  MASS  AND  DAMPING  COEFFICIENTS  FOR  u  -  o. 


As  the  frequency  of  the  motion  b>  approaches  to  zero  the  values  of 
m^  given  by  (2.22)  goes  to  zero  as  well  while  m^h  values  tend  to  qir. 

The  coefficients  of  the  linear  equation  (2.36)  become  real  and  the 
formulation  correspond  to  the  case  where  rigid  boundary  conditions  exist 
at  the  free  surface.  The  only  imaginary  value  Is  the  imaginary  part  of 


Aq  which  Is  equal  to 


Im  A^  =  A.^  = 


(4.1) 


From  which  using  3.9  the  damping  coefficient  Is  immediately  obtained  as: 


11m  “22 

u)  0  (ijpD 


.a2 

4h(h-d) 


(4.2) 


The  above  result  can  be  also  obtained  by  using  the  Haskinds  relations 
given  for  finite  depth  by  Newmann  (1976).  The  equat1on(2.36),(2.37)and 
^.38) for  heave  motion  are 


,0  ■  "o<V> 


0°  E 
^0  00 


q'l 


«  H-(m-a) 

"  Ho'(m^a) 


-  (y  ’  1 


(4.6) 

(4.7) 


The  above  equations  for  very  small  m^  values  can  be  written  as  follows 
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(4.8) 


Yoj  =  <Ji(0,j)  +  i  e  (0,j) 

^nj  “  ^ 


(4.9) 


where 

''J  ^  T  i  7/  I''  -  IL.  \ 


q*l  n 


q  (sin  (qird/h))^ 


lir  K  -  nr 


(4.10) 


o  ah^  ?  K  (qua/h)  ,  ?  i  i  ,  2  ,2 

’o  =  ^  Kf(qTa/h)  ‘  ^  qnd/h)^-(]4(l)  >>  “  ^ 


h  .  {Zihi,.,,"  f  .  I  (sii 

"n  ,3j  CTq.a/h)  q  [tq3 

q=i 


1  (sin  qird/h] 


.  2MA  ^ 

(nTr)2  ^  ^  2it 


iDU^oiy 


(4.11) 


After  these  preparations  taking  the  limit  of  (4.4)  for  very  small  m^^ 
values  and  using  (4.8),  (4.9),  (4.11)  we  obtain 

J  (V  ^  ^i.)  "’o  ®  =  h^^  +  i  m^  h.^  (4.12) 

J  *  1  J  j 

Here  subscript  r  and  i  refer  to  the  real  and  imaginary  part  of  the  values. 
Separating  into  real  and  imaginary  parts 

?  n>J  3  (n^j)  +  A.j  a  (n,j)]  =  h.^  m^ 

z  [A^j  ot(n.J)-A^.j  mj  $(n,j)]  =  h^^. 

Taking  the  limit  of  the  right  hand  side  of  the  first  equations  as  m^-*-  o 
we  obtain 

.f  [a  (n,j)  +  mj  B  (nj)]  =  h^^ 

J 

and  A^j  a(n  J  )  *  -  mj  6  (n,j  )  A^j 
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One  can  see  therefore  that  all  A.'  s,  j  *  1,2...  have  real  values. 

As  a  special  case  let's  take  a  floating  disc.  This  can  be  expressed 
as  a  limiting  case  where  d^h.  One  can  show  that  for  this  case  the  equation 


4.4  has  a  diagonal  matrix  and  the  unknown  A^'  s  are  given  by 


1  /a. 
7 


(4.13) 


Since  all  sine  terms  In  equation  4.12  are  equal  to  zero.  The  added  mass 

of  tht  circular  disc  In  heaving  motion  for  a)-»-o  Is  obtained  as 

11m  **22  _  /hx  ,  1  .  1  /a,  4  /h.  «  1  ^1  (4.14) 

^  (?>  ^  5  <ir)  -  P  (?) 


In  another  limiting  case  where  d  goes  to  zero,  again  all  the  sine 
terms  In  equation  4.12  are  equal  to  zero  and  the  coefficient  matrix  becomes 
diagonal.  In  this  case  one  can  show  that  as  u  and  d  go  to  zero  the  heavi 
added  mass  In  fact  Is  Infinite.  According  to  this  present  theory  one  can 
conjecture  that  for  d<h  and  at  finite  h  the  heave  added  mass  remain  finite 
at  low  frequencies  but  qoes  to  Infinity  for  Infinite  depth. 


5,  NUMERIC  SOLUTION  AND  RESULTS. 

The  specific  formulas  used  for  the  calculation  of  added  mass  and  damping 
coefficients  are  given  In  the  Appendix.  The  required  routines  for  the 
calculation  of  Bessel  functions  and  solution  of  linear  equation  are  obtained 
from  IMSL  computer  library.  The  computations  were  done  at  the  U.  S.  Naval 
Academy*  The  results  for  heave  were  first  compared  to  the  results  published 
by  Garrison  (1975)  and  Kritls  (1979).  The  results  are  given  In  Figure  2. 
Added  mass  values  obtained  by  this  theory  compared  well  with  those  of 
Kritls  while  Garrison's  number  are  observed  to  be  higher.  The  damping 
coefficients  for  heave  are  observed  to  be  less  than  the  values  reported 
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by  Kritis  (1979)  and  they  were  observed  to  agrte  better  with  those  of  Garrison. 

The  heave  added  mass  and  damping  coefficients  are  also  compared  with 
the  experimental  data  reported  by  McCormick,  et  al  (1980). 

Figure  3  shows  experimental  data  and  theoretical  values  computed  for 
0)  =  3  rad/sec.  In  this  diagram  added  mass  value  is  nondimensionalized  by 
the  total  mass  of  the  cylinder  of  height  equal  to  the  depth  of  water. 

This  is  expressed  in  the  diagram  as  ANU/MH.  The  experimental  data  are 
observed  to  remain  above  the  theoretical  curve  the  best  correlation  is 
observed  at  about  depth/radius  values  equal  to  7. 

Experimental  and  theoretical  damping  coefficients  are  compared 
in  Figure  4.  Experimental  values  remained  above  the  theoretical  calculation 
while  showing  similar  trends.  This  discrepancy  can  be  possibly  explained 
by  the  viscous  damping  neglected  by  the  theory.  The  numerical  results  are 
also  tested  with  those  reported  by  J.  Bai  (19/6)  and  the  results  for  added 
mass  and  damping  coefficients  are  observed  to  agree  with  in  a  derivation  of 
3  percent.  In  all  calculations  the  infinite  series  are  represented  by 
20  terms. 

The  heave  added  mass  values  for  different  water  depths  and  draft  values 
are  given  in  Figures  5  to  HI.  The  general  behavior  of  the  curves  is  that  as  the 
frequency  increases  the  values  remain  constant.  At  zero  frequency  numerically 
at  least  the  added  mass  values  are  observed  to  increase.  The  behavior  of 
the  curves  at  small  frequency  is  discussed  in  section  4.  Damping  coefficients 
for  heave  are  given  in  Figures  12  to  17.  At  high  frequencies  these  values 
are  observed  to  tend  to  zero  while  at  zero  frequency  the  values  are  finite 
at  shallow  depth  and  tend  to  approach  zero  as  the  depth  increases.  Higher 
values  are  observed  to  correspond  to  small  draft  to  radius  values.  Deep 
water  cases  correspond  to  h/a  =  20. 
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The  sway  calculations  are  compared  to  the  results  published  by 
Bai,  Yeung  (1974)  and  are  presented  in  Figure  19.  This  computation  tends  to 
follow  the  values  computed  by  Bai  while  the  values  reported  by  Isshiki 
remained  low  especially  at  peak  values.  Figures  20  to  25  ^how  the  sway 
added  mass  values.  Added  mass  values  are  observed  to  increase  as  draft 
to  radius  is  increased  at  all  depths.  These  values  also  remained  finite 
at  small  frequencies.  All  curves  are  observed  to  have  a  local  maximum  at 
about  V  =  1.  Damping  coefficients  for  sway  are  equal  to  zero  at  zero 
frequencies  and  are  observed  to  increase  as  the  draft  increases.  The  curves 
for  damping  coefficients  are  given  in  Figures  26  to  31.  At  high  frequencies 
the  curves  show  a  decreasing  slope  and  the  maximum  values  are  again  observed 
at  about  v  =  1.  Pitch  computations  are  first  checked  with  those  reported  by  Bai 
and  Garrison  and  a  good  agreement  is  observed. 

Pitch  moment  of  inertias  are  presented  in  Figures  32  to  37.  Except 
at  very  small  drafts  the  curves  have  a  very  small  slope.  Inertia  coefficients 
are  observed  to  decrease  as  draft  increases  at  shallow  waters  while  in 
deep  water  higher  coefficients  correspond  to  higher  drafts.  Figures  38  to 
43  give  thepitch  damping  coefficients.  At  shallow  water  (h/a  =  1  h/a  *  3) 
high  damping  coefficients  correspond  to  small  drafts  while  at  deep  water 
(  h/a  =  20)  high  damping  coefficients  are  observed  to  correspond  to  high 
drafts.  A  peculiar  curve  is  seen  in  Figure  42  for  T/A  =  0.1  which  suggests 
that  the  damping  coefficient  for  very  shallow  disks  increases  as  draft 
decreases  even  at  moderately  deep  waters.  Pitch  induced  sway  added  mass 
(a^i)  and  pitch  induced  sway  damping  coefficients  bgi  are  given  fn  igure  44-46. 
It  is  interesting  to  note  that  some  of  these  values  are  in  fact  negative 
for  low  drafts  at  finite  depths  such  as  h/a  *3  but  as  the  draft  increases 
the  values  become  positive. 
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CONCLUSION 


The  solution  presented  in  this  paper  offers  a  quick  calculation  of 
the  hydrodynamic  coefficients  for  a  simple  verticW. circular  cylinder. 

The  results  are  compared  to  some  available  experimental  numerical  results. 
The  agreement  with  numerical  results  are  observed  to  be  satisfactory.  The 
comparisons  with  the  experimental  showed  that  even  though  the  trend  is 
well  represented  the  amplitudes  are  not.  This  is  partially  acceptable  at 
least  for  the  case  of  damping  coefficients  where  the  viscous  resistance 
must  be  effective.  The  input  variables  are  water  depth,  radius  of  cylinder 
and  its  draft.  Well  documented  routines  can  be  used  for  the  calculation  of 
special  functions  and  the  solution  of  the  linear  equations.  The  necessary 
formulas  for  computer  application  are  also  presented. 

The  hydrodynamic  coefficient  to  study  the  motion  of  the  cylinder  are 
presented  in  graphical  format.  It  is  hoped  that  this  will  increase  the 
efficiency  of  future  designs  and  that  the  designer  will  be  able  to  estimate 
these  coefficients  rather  precisely  for  his  computations. 

The  limit  of  the  heave  added  mass  value  for  zero  frequency  is  also 
discussed.  It  is  shown  that  this  quantity  remains  finite  for  finite  depth 
and  goes  to  infinity  for  infinite  water  depth.  Special  formulations  are 
seen  to  be  required  for  this  limiting  case. 

The  present  formulation  is  currently  extended  to  cylinders  with  vari¬ 
able  cross  sections. 


APPEHDIX 


FORMULAS  FOR  COMPUTER  APPLICATIONS 


A>Heave 

The  linear  set  of  equations  for  complex  coefficients 

Y  nj  ''j  “  K 

is  solved  first.  The  expression  for  and  h^^  are  given  below 


■  (V’‘  “j*"-  "o'"-”*',  Pq'"-’”  *  *  “ 


32  UjTPp(nj) 


m, 


oa 


00  Ko(mna)  r  ft 


wliere 


U 


,  =  ^  j  =  1.  2.  3  .  . 


(-l)"'^'^  (m  d  sh  (m  d))^ 

P-(n»j)  ~  . .  A  "  n  9  1 1  ^ 

°  [(2m  h)  +  sh(2m  h)lt(m  d)^  +  (rrr  )‘^]  [(m^d)^  *{5^r] 


(-1)"'^^  (m^d  sin(mqd))^ 

[(2m  h)  +  sin(2m  h)][(m„d)^  -  (nn)^]  [(m„d)^  - 


Jo(moa)  (m^a)  ^  Yo(moa)  Y^(mpa) 

[Jl(moa)]2  + 

_  1 _ 

[0l(moa)]2  +  [Y,(m^,a)J^ 
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The  added  mass  and  damping  coefficients  are  then  computed  by  the  following 
equation: 


!22  +  !22  ,  _cL  ri.  i ^2  ^  A  (-1)"  1^(5^) 

pD  ^  O.PD  h-d  ^2  8  ^d)  "  2  ''o  ^ V-  r;^  > 

where  ^22’  ^22  added  mass  and  damping  for  heave,  D  is  the  volume 

of  the  cylinder  in  water. 


B.  Sway 

The  linear  set  of  equations  for  sway  can  be  written  in  the  form: 

I'nJ  ■  -  ''n 

where 

..  32/ 


.  4a-P„(n,o)K„  D,-P,(n.o)K,)  -  i^T.p„(n.o)K„ 


and  where 
U- 


I  (j  ^a\ 

'j  '  ■  A’  “o  '  fi 


2  2,22 


T-  (CJ„("i„ii)  *  V„(v)l  ♦  liV  '■'I'V)  *  ’'ifv"  -  ST 

0  0 

NotV^'^l^V^  +  YQ(mQa)Y^(mQa)]f^ 


1  ^  2 

L  =  {[Jo(mQa)J^(m^ja)  +  Yg(mQa)  Y^Cm^a)]  -  —  [J^Cm^a)  +  Y^^Cm^ja)]}! 

0 

sh(m-h)  -  sh(m  d)  sin(m  h)  -  sin(m  d) 

<0  •  (  -WdV  ■  >  •  S  •  ' - stMipj')  ^  > 


“k  ■  - 


m  a 

^  K  fm  ai  ;  ^  ^ »  2,  3,  ,  .  . 

1  +  O^V^ 
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The  definition  for  and  remain  the  same  as  In  heave. 

The  sway  added  mass  and  damping  coefficients  are  obtained  as 


n  Ish(m  h)»sh(m  d)]  »  K,(m,^a)(s1n(m  h)-s1n(m  d)] 

1  -  {B  H,(m  a)  5 - 5 — +2:  B„  '  — 3 - 9— 

D  (i>pD  0  1  0  N  ^  m  (h~d)  ^  u  ^  —  /i.  j\ 


pD  ujpD 

where 


q*l 


mq(h-d) 


®q  '  hmqFq(a)  ^4a  ^qo  ^n  "^n  ^qn^  *  hm^  F"(a) 
n  *  0,  1,  2,  ..  q  =  0,  1,  2,  ... 


_  (-1)\~^mQd  sh(mpd) 
[(m.d)^  +  (nii)^] 


■qn 


(-1)”Nq"^2mqd  sinCm^d) 

((m^d)  -  (nir)^] 


In.-**  In:**  [Sin(m  h)  -  sin(m  d) 

0  a  0  ‘ - - -  ’  pq  a  q  - ^ ^ — ■■ 


and 


»  1  /I  .  s1nh(^moh  ) 


,  s1n(2m  h) 

Nq  “  ?  (1  + - ^  ) 

3  2  2m^h 
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C.  PITCH 

The  linear  set  of  complexed  valued  equations  is 


^nj  %  n  =  0,1.2  ,  j  =  0,  1,  2  ... 


where 


T„  ■  (  -  5„j  *  U  Uj  [L-P„(„.j)  0,-P,(r,J)))  -  1  ^UyT.P^(„.j) 


\  =  8  (L-P„(n,o)-  t  E  0,-P(|(n.o)-e|,>  -A  -1  ^  -T-P^Cn.o)  ^ 

q=l  ^  ^  0 

I  (  ra/d)  ,  ft  • 

“J  =  T^TTiTdT'  “o  '  ^  •  ‘nj  *  <1  V-l 

_ 3  h  3,a,2  .  I  .  .  cn(m^h)-2ch(m„d) 

®o "  ■  d  •  8^d^  ^  ^  ^ 

_ 3  h  3.ax2  1  -  2ch(mqd) 

®q  "  -  TS^l  -  {i^dTT-ln  {h^Sy  ^ 


>forn  =  1.2,  3.  ... 


“o^“  ^  1  ■  I  }  for  n  =  0 


The  expressions  for  P^,  P^,  L,  T,  remain  the  same  as  defined  for  sway. 
The  roll  (Pitch)  added  moment  of  inertia  is  calculated  as 


’66  ■"  '  ‘j6/^ - d_  ,  i  H,  (n.„r)  B„  +  E  t„  K,  (m„a)  B„  J  + 


pb  a2 


h-d 


r  (1-1  (i)  )  +  l/i)  A  +  1  E  ^'1^  — _ ^nl> 

I  6  ^d^  ^  Ird^  %  IT  n  I^(nra) 
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The  expression  for  B^,  remain  the  same  except  the  B*  values  for  roll  are 


given  by 


3  h  3  »  2  1  sh(m  d)  ch(ni  h)  -  2ch(«i  d) 


o  h  a  1  \—n"/  “  ~ — '"'n' 

“  (m-d)iJ  ^  (m'3)  ^(m  d)2  ^  ^ 


sin(mpd)  cos(m„h)-2  cos(m„d) 

V 


and  the  expressions  for  t^  and  t^  are 


to  =  ^ 


^  sh(mQd)  ch{mQh)  -  ch(mjjd) 

'  — rar? - ^ 


tq  =  -  ?)  ^ 


0  '0 
s1n(m  d)  cos(m„h)  -  cos(m  d) 
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